Dual quark condensate in the Polyakov-loop extended Nambu-Jona-Lasinio model 
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The dual quark condensate E^"^ proposed recently as a new order parameter of the spontaneous breaking 
of the Z3 symmetry are evaluated by the Polyakov-loop extended Nambu-Jona-Lasinio model, where n are 
winding numbers. The Polyakov-loop extended Nambu-Jona-Lasinio model well reproduces lattice QCD data 
on X''-^^ measured very lately. The dual quark condensate S^"^ at higher temperatures is sensitive to the strength 
of the vector-type four-quark interaction in the Polyakov-loop extended Nambu-Jona-Lasinio model and hence 
a good quantity to determine the strength. 
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Recently, a new order parameter of the Z3 center sym- 
metry was proposed by using the chiral condensate a and 
evaluated by quenched |1, 2] and full lattice QCD J3t]; the 
new order parameter is called the dual quark condensate. 
This makes it possible to discuss the connection between 
the quark confinement and the chiral symmetry. The rela- 
tion can be discussed also by the Polyakov-loop extended 
NJL (PNJL) model |4] in which the Polyakov loop <P is ap- 
proximately treated as a classical variable. Actually, ex- 
tensive studies are made on the relation; for example, see 
Refs. illlSI?, 8, 9, 10,[ni[il[il[il[lll- The dual quark 
condensate is calculated also in the color SU (2) system by the 
Dyson-Schwinger equation |16], but this approach does not 
treat the confinement mechanism dynamically. Very recently, 
the functional renormalization-group method was also applied 
to evaluate the dual quark condensate in the color SU (3) sys- 
tem lO. 

In this paper, we evaluate the dual quark condensate by us- 
ing the PNJL model and show that the PNJL result can re- 
produce full lattice QCD (LQCD) data (13] on the dual quark 
condensate. 

We consider the quark field q that obeys a twisted temporal 
boundary condition 



g(x,/?) 



=g(x,0), 



(1) 



where (p is a twisted angle. Now we define the (/3-dependent 
chiral condensate Illl|j,l3|] by 



f7M = --(Tr[(,n + D^)-i]), 



(2) 



where the twisted boundary condition ([T]i is imposed on the 
Dirac operator D^p, but the bracket (■■■) keeps the anti- 
periodic boundary condition (/? = tt IH [31] . The dual quark 
condensate S'^"^ is defined by 



Ztt 



(3) 
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with winding numbers n d Silt]. In particular, Z'^^^ is called 
the dressed Polyakov loop. The winding number is 1 in both 
Z'(i) and In full LQCD calculations of Ref 13], the dressed 
Polyakov-loop is evaluated with m/T = 0.032. We can ex- 
pect that the two have similar T dependence to each other. 
Thus, the dual quark condensate relates the chiral condensate 
to the quark confinement iflUlUfTelfTTifTsI]. 

We start with the two-flavor PNJL Lagrangian with the 
vector-type four-quark and the scalar-type eight-quark inter- 
actions; see Ref. ]8] for the details. The scalar-type eight- 
quark interaction makes the chiral phase transition stronger 
Hence, the PNJL result becomes consistent with LQCD data 
at finite T ]& JJ2, J9J. Furthermore, both the scalar-type eight- 
quark and the vector-type four-quark interaction are necessary 
for the PNJL model to reproduce LQCD data ll20l[2lil22|] at 
imaginary chemical potential [lOj . 

Making the mean field approximation to the PNJL La- 
grangian and performing the path integration of the resultant 
partition function over q under the twisted boundary condi- 
tion ([T]i, one can obtain the thermodynamic al potential 



n = -2Nf 



1 



(2^ 



3E„ 



^ In [1 + 3{<P + <ie-''^p )e-''^p 



1 



+ - In [1 + 3(<Z> + <Pe-'^^^ ) 



-3/3iJ- 



(4) 



where a = (gg), uj = (q^oq) and 
with Ep = v/p2 + M2, M = Too - 2Gsa - 4Gs8cr3 and 
iT9^ = —2GvOJ — iT^T + iTip. Here, toq is the current quark 
mass and we use toq = 5.5 MeV. The constants Gg, Gv, Ggs 
denote coupling strengths of the scalar-type four-quark, the 
vector-type four-quark and the scalar-type eight-quark inter- 
action, respectively. The Polyakov-loop potential is a func- 
tion of and its Hermitian conjugate ^. The 3-dimensional 
momentum integration is regularized by a cutoff A. The clas- 
sical variables X = a,uj,(I>, <P are determined by solving the 
stationary conditions dfl/dX ~ numerically, and Z'*^"' is 
obtained numerically from (7{^) with ([3]). Details of the PNJL 
model are shown in Ref. liSJ . 



2 



In the PNJL calculations, two types ofU are often used; one 
has the polynomial form |5] and the other has the logarithm 
form [6]. The former (latter) potential is referred to as RTW05 
(RRW06). Parameters of these potentials ai-e fitted to LQCD 
data at finite T in the pure gauge limit |23, 24], but one of the 
parameters, Tq, is usually refitted to reproduce the pseudo- 
critical temperature Tc = 173 ± 8 MeV in full LQCD Q- 
The value of Tq thus determined is 185 (200 MeV) for the 
RTW05 (RRWOeVrntential. The RRW06 potential can repro- 
duce LQCD data i20ll2ll l22^ at imaginary chemical potential, 
but the RTW05 does not [10]. Therefore, the RRW06 poten- 
tial is mainly used in the present PNJL analyses. 

We do four types of PNJL calculations shown in TableU In 
PNJL-I, Gs and A are adjusted to the pion mass Af^ — 138 
MeV and the pion decay constant and f„ = 93.3 MeV at 
T = 0; see Ref. |t8J for values of the parameters. In PNJL- 
II, -III and -IV, Gs, G^s, A are fitted to AU = 138 MeV 
and = 93.3 MeV and the sigma meson mass M^r = 600 
MeV; see Ref. [8] for values of the parameters. Since the 
empirical value of M^r has a large error bar, the empirical 
value is not a strong constraint to determine the strength Gss- 
We then think Gss as a nearly-free parameter The value of 
Gs8 determined above is just an example. In PNJL-III, the 
vector-type interaction is a free parameter The strength of the 
vecto r-type i nteraction was determined to reproduce LQCD 
data 1I20I I211 12?! at imaginary chemical potential |10]. The 
ratio Gv/Gs is 0.85. As an example, we take a bit smaller 
value of Gv/Gs = 0.667. For comparison, we also do the 
NJL calculation with the scalar-type eight-quark interaction; 
this is referred to as NJL. Since NJL agrees with PNJL-II in 
the limit of T = iQ], the same parameter set is taken in the 
two models. 



Model Interaction U rf'" line 

PNJL-I RRW06 LOlTc 1.32rc l.OlTc thin-solid 

PNJL-II (7^ RRW06 Tc 1.14rc Tc solid 

PNJL-III a^ a^, cj^ RRW06 Tc 1.14Tc Tc dashed 

PNJL-IV (T^ RTW05 0.98Tc l.lOTc l.UTc dotted 

NJL o-^ a* - - 0.94rc 0.97rc dot-dashed 



TABLE I: Definitions and results of PNJL and NJL calculations. 
Symbols, RTW05 and RRW06, are the Polyakov potentials of poly- 
nomial type |5] and logarithm type |6], respectively, while a"^, cr*, uj'^ 
denote the scalar-type four-quark, the scalar-type eight-quark and the 
vector-type four-quark interaction, respectively. Symbols T" , T*, 
denote the pseudo-critical temperatures defined by peak posi- 
tions of da/dT, d<P/dT, dr*^' /dT, respectively. Here, Tc = 173 
MeV. 



As mentioned above, the twisted boundary condition ([T]i is 
imposed on D^p . Therefore, <P and <P are first obtained under 
the anti-periodic boundary condition (p — tt. The quantities 
(j{(p) and uj{ip) are determined by solving the stationary con- 
ditions, dn/d(j — and dfi/du) = 0, numerically under 
the twisted boundary condition ([T]i with and ^ fixed to the 
values determined alip — ir. 

The pseudo-critical temperature Tf of the deconfinement 



crossover is usually defined by a peak position of either the 
Polyakov-loop susceptibility or <M?/dT. In the present analy- 
sis, we take the latter to compare Tf with the pseudo-critical 
temperature T^^ ' estimated by a peak position of dS^^^ /dT, 
because we can not define the susceptibility for the dressed 
Polyakov loop For consistency, the pseudo-critical tem- 
perature of the chiral crossover is defined by a peak po- 
sition of da/dT. In the PNJL calculations with the RRW06 
potential, da/dT has two peaks, but we take the second peak 
close to a peak position of the chiral susceptibility. 

TableUsummarizes values of T*, T^, T^^ in five types of 
model calculations, where the values are normalized by Tc — 
173 MeV, i.e. the LQCD result f23\ for Tf and T^ . In PNJL- 

I, -II, -III with the RRW06 potential, Tf and are close 

to each other and also to the LQCD data ll25ll . although this 
property is not seen in PNJL-IV with the RTW05 potential. 
As for Tc, PNJL-II is more consistent with the LQCD data 
than PNJL-I and hence the scalar-type eight-quark interaction 
is necessary to reproduce the LQCD data. 
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Fig. 1: Chiral condensate as a function of T /T^ and renormalized 
Polyakov loop as a function of T/Tf in the case of ip = TT. Here, 
o is normalized by the value ao at T = 0. The values of and 
T* are summarized in Table. U Definitions of lines are also shown 
in Table. H] LQCD data shown by box and triangle symbols are taken 
from Ref. and Ref (23], respectively. 

Figure [U presents T-dependence of the chiral condensate 
and that of renormalized ^, where the anti-periodic boundary 
condition = tt is taken. The resultant chiral condensate 
cr is normalized by the value cro at T = 0. Temperature T 
is normalized by T^ for the chiral condensate and by Tf for 
renormalized to compare the shapes of a/a^ and <P with 
LQCD data Also for the shape for ct, PNJL-II (solid 

curve) gives a better agreement with the LQCD data than 
PNJL-I (thin solid curve). Meanwhile, the solid and thin solid 
curves are almost identical for 'P, indicating that the scalar- 
type eight-quark interaction hardly affects the shape of 'P. At 
zero quark chemical potential, the vector-type four-quark in- 
teraction does not affect a and (p in the mean field approxi- 
mation; hence, PNJL-III yields the same a and (p as PNJL-II 
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in Fig. [U Nevertheless, the interaction gives a sizable effect 
on as shown later For the shape of <P, PNJL-IV (dotted 
curve) gives a larger disagreement with the LQCD data than 
PNJL-II. Thus, the RRW06 potential is better than the RTW05 
potential also for the shape of (P. 

Full LQCD simulation of Ref. H shows that T^^*" ~ 
~ Tf ~ 153 MeV, while that of Ref. (lE] does T^^ ~ 
T* ~ 173 MeV. The two LQCD results have a systematic er- 
ror of ^ 10 %. Moreover, it is reported in Ref. lEsll for the 
2 + 1 -flavor system that there exists a non-negligible devia- 
tion between and T*. This indicates that it is an unsettled 
problem whether and T* really coincide or not. In this 
work, however, we assume the coincidence, since the LQCD 
data 1 3] that we are analyzing has the property. In the present 
analysis, temperature is normalized by T^, but it is taken to be 
153 MeV for LQCD data of Ref. H and 173 MeV for model 
calculations and LQCD data of Ref. ll25l[26il27ll. 

Figure |2] presents the normalized dressed Polyakov loop 



^(^V^T? as a function of T/Tf 



where s!^'' is the 



dressed Polyakov loop at T T^-^ . PNJL-I (thin solid line) 
and PNJL-II (soUd Une) well reproduces LQCD data 13] (box 
symbols), but PNJL-IV (dotted line) and NJL (dot-dashed 
line) do not. Near T^^^\ the scalar-type eight-quark inter- 
action hardly affects the dressed Polyakov loop. Thus, the 
RRW06 potential is necessary to explain T dependence of (p, 
a and S'^^^ consistently. 
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Fig. 2: Normalized dressed Polyakov loop as a function of r/r<f . 
Definitions of lines are shown in Table. |l] LQCD data (box symbols) 
are taken from Ref. [30- 

Figure [3] shows (^-dependence of cr at T = 200 and 250 
MeV, while Fig. |4] represented T-dependence of the dual 
quark condensates with rt = and 1. The difference between 
PNJL-I (thin-solid line) and PNJL-II (solid line) shows an ef- 
fect of the scalar-type eight-quark interaction. As shown in 
Fig.[3j the effect on a is large around Lp = tt, but small around 
tf = Q. The absolute value of a is relatively larger around 
tf = Q. As a consequence, as shown in Fig. |4l the effect on 



0.04: 



> 



"D 0.02 





1 =£iou [ivievj 














-T=200 




[MeV]\>^;x 
* \^\\ 


/ /m 1 


\ \\\\ 


///I 1 




y /' 1 1 




// I 










/' / / / 


\ \ V 

\ \ \ 


//'' 


\ \ \ 


/ / / 


\ \ V 


y / / 


\\ 


// 







1 

(p/7[ 



Fig. 3: (/3-dependence of the chiral condensate at T 
250 MeV. Definitions of lines are shown in Table. U] 



0.03' 



200 and 



S 0.02 



0.01 




0.2 0.4 
T [GeV] 



0.6 



Fig. 4: T-dependence of the dual quark condensate. Definitions of 
lines are shown in Table. |l] 



and are appreciable only at T/^*" < T < 2T^''^\ 
Thus, the effect of the eight-quark interaction becomes more 
appreciable for tr at = tt than for 17'^°) and Z"'^'. At 
T > 2T^ ^ , the effect becomes negligible, since a itself 
is tiny there. Meanwhile, the difference between PNJL-II 
(solid line) and PNJL-III (dashed line) presents an effect of 
the vector-type four-quark interaction. The effect on a is zero 
at = TT, and appreciable at nonzero tt), as shown in 
Fig. [3] The effect on and Z'^^^ is then appreciable at 



Ti^*" < T < 1.5Ti^*" and sizable at T > l.bT^'^\ 



as 



shown in Fig.|4] Thus, the effect of the vector-type four-quark 
interaction becomes more appreciable for S'^^^ and Z^^^ than 
for cr at </j = TT. This indicates that the dual quark condensate 
at r > 1.5T(,'^* ' is a good quantity to determine the strength 
of the vector-type four-quark interaction. 

The present PNJL results are consistent with the LQCD 
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data for 17^^^ / near Tc, but not for (j{ip) itself even if 

both are compared at the same values of T/Tf'^ \ although 
the latter is not shown explicitly in this paper. In LQCD, (T{(p) 
and 17*^^^ are not renormalized, while presented in Fig. [T] 
is renormalized. The reasonable agreement of the PNJL result 
with the LQCD one for ct/ctq and <P in Fig.[I]and S'^'^') / S^^^ in 

Fig . |2] may imply that cr/ ctq and 17^^^/ ^ have better renor- 
malization properties than ij{ip) itself. 

Finally, we discuss the sensitivity of S^^^ to the parameters 
of the PNJL model. If A is varied by 10 %, the pion decay 
constant at T = is changed by about 20 %. Hence, one can 



change A by only 0. 1 % in order that the calculated pion decay 
constant is consistent with the observed one with 0.1% error 
This is also the case for Gg. If and Gg are varied by 1 %, 
are changed by about 10 %, but the normalized quantity 
S^^^ 1 17^^^ is hardly changed near Tc. Further, Z'(i) near is 
much less sensitive to Ggg and Gv, although these are nearly- 
free parameters. Thus, we can think that T,'^^^ / near Tc 
changes little from the present values. 
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